In this case study, we will build an integrated scRNA-seq analysis workflow to interpret the transcriptional and clonal structure of tumor-infiltrating T cells in MC38 colon adenocarcinoma (data from Xiong et al 2019).

The main R packages and methods employed in this workflow are:

Note that scRepertoire requires R version 4.0 or higher - you will need to have R v.4 installed to run this case study.

scTCR data preparation

Download scTCR data from Array Express (https://www.ebi.ac.uk/arrayexpress/experiments/E-MTAB-7918/)

Mouse 1 to 4 (sample ID 4 to 7) correspond to TCR-seq libraries 35 to 38

Combine alpha and beta chains using the combineTCR function from scRepertoire

The function ‘combineExpression’ of scRepertoire allows incorporating clonotype information to a Seurat object, and creates the categorical variable ‘cloneTypes’ discretizing frequencies of the clonotypes

We have now paired expression and TCR data for the query samples, and loaded them into a unified Seurat object. We can proceed to project the data onto the reference atlas.

ProjecTILs

Load the reference atlas

[1] "Loading Default Reference Atlas..."
[1] "/Users/mass/Documents/Projects/Github/ProjecTILs_CaseStudies/ref_TILAtlas_mouse_v1.rds"
[1] "Loaded Reference map ref_TILAtlas_mouse_v1"

Project query data (with clonotype information stored as metadata) onto the TIL reference atlas

[1] "Using assay RNA for query"
celltype.pred
         CD4T          CD8T     Non-Tcell Tcell_unknown          Treg 
          467          3115           431           322          1739 
      unknown 
          108 
CD8Tstate
         Naive EffectorMemory     MemoryLike      Exhausted        unknown 
           180            754            570           1510            101 
[1] "539 out of 6182 ( 9 % ) non-pure T cells removed.  Use filter.cells=FALSE to avoid pre-filtering (NOT RECOMMENDED)"
[1] "Aligning query to reference map for batch-correction..."

Projecting corrected query onto Reference PCA space
[1] "Projecting corrected query onto Reference UMAP space"

Visualization of projected data.

T cells are projected mostly in the Treg, CD8 Terminal exhausted (CD8_Tex) and Precursor exhausted (CD8_Tpex) areas, with some clusters in the T-helper (Th1) and CD8 Effector memory areas, and to a lesser extent in Naive-like and CD8 Early-Activated areas.

Visualize the projections per sample. Broadly, the distribution across the atlas is similar for the four mice, with some variation in the frequency of Effector Memory T cells.

Classify projected T cells into cell subtypes/states


     CD4_NaiveLike     CD8_EarlyActiv CD8_EffectorMemory      CD8_NaiveLike 
                31                109                575                261 
           CD8_Tex           CD8_Tpex                Tfh                Th1 
              1752                436                 10                541 
              Treg 
              1928 

Look at distribution of T cells in terms of cell states.

We can check the gene expression profile of the cells assigned to each state (yellow), and compare them to those of the reference states (black).

For example, cells projected into Tex and Tpex region express Tox and Pdcd1, while only Tex express Gzmb and Havcr2, as expected.

Clonality analysis

Now, let’s see where the expanded clones are located within the map:

Larger (most highly expanded) clones are concentrated in the Tex/Tpex area. This is largely explained by the fact that sustained antigenic stimulation drives immunodominant clones into the Tox+-driven exhausted lineage (Tex - Tpex region)

Let’s check the raw numbers

                                
                                 CD4_NaiveLike CD8_EarlyActiv
  Hyperexpanded (100 < X <= 500)             0              0
  Large (20 < X <= 100)                      0              1
  Medium (5 < X <= 20)                       0              2
  Small (1 < X <= 5)                         0              5
  Single (0 < X <= 1)                       23             42
                                
                                 CD8_EffectorMemory CD8_NaiveLike CD8_Tex
  Hyperexpanded (100 < X <= 500)                  1             0      86
  Large (20 < X <= 100)                          99             0     405
  Medium (5 < X <= 20)                           47             1     313
  Small (1 < X <= 5)                             91             5     216
  Single (0 < X <= 1)                           129           123     138
                                
                                 CD8_Tpex Tfh Th1 Treg
  Hyperexpanded (100 < X <= 500)       31   0   0    3
  Large (20 < X <= 100)                98   0   0    8
  Medium (5 < X <= 20)                 68   1  27  121
  Small (1 < X <= 5)                   54   2 118  225
  Single (0 < X <= 1)                  32   3 170  287

Plot clonal size by functional cluster

This confirms the intial observation that Tex and Tpex present the largest clonal expansion

We can highlight specific clonotypes on the reference atlas (here those with at least 40 cells):

The majority of clones tend to span Tex and Tpex states. Indeed, Tcf7+ precursor exhausted Tpex cells self-renew and give rise to more differentiated Tex effector cells.

The clonal overlap (Morisita similarity index) implemented in scRepertoire confirms this pattern:

We can also visualize overlap of the top CD8 clones, and their overlap across functional clusters, using scRepertoire’s alluvial plot

Conclusions

Projection of the MC38 TILs onto the reference atlas showed that these T cell samples mostly consist of exhausted (CD8_Tex), precursor-exhausted CD8 T cells (CD8_Tpex), and Tregs, typical of highly immunogenic, “hot” tumors. Also, the majority of expanded clones were found spanning the CD8_Tex and CD8_Tpex states. This is expected as sustained antigen stimulation in the tumor drives immuno-dominant tumor-reactive clones towards the exhaustion differentiation path.
The combination of ProjecTILs and scRepertoire simplifies the joint analysis of single-cell expression data and clonotype analysis, in the context of an annotated reference atlas of TIL states.

Further reading

Original publication - Xiong et al. (2019) Cancer Immunol Res

ProjecTILs case studies - INDEX - Repository